physica status solidi, 23 September 2009 



Reentrant phenomenon in the 
exactly solvable mixed spin-1/2 
and spin-1 Ising-Heisenberg model 
on diamond-like decorated planar 
lattices 

Lucia Canova 1 '*, Jozef Strecka 2 

1 Department of Applied Mathematics, Faculty of Mechanical Engineering, 
Technical University, Letna 9, 042 00 Kosice, Slovak Republic 

2 Department of Theoretical Physics and Astrophysics, Faculty of Science, 
P. J. Safarik University, Park Angelinum 9, 040 01 Kosice, Slovak Republic 

Received XXXX, revised XXXX, accepted XXXX 
Published online XXXX 

PACS 05.50.+q, 75.10.Hk, 75.10.Jm, 68.35.Rh 

* Corresponding author: e-mail lucia.canova@tuke.sk, Phone: +421 55 602 2228, Fax: +421 55 633 4738 



Ground-state and finite-temperature behaviour of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on decorated 
planar lattices consisting of inter-connected diamonds is investigated by means of the generalised decoration-iteration 
mapping transformation. The obtained exact results clearly point out that this model has a rather complex ground 
state composed of two unusual quantum phases, which is valid regardless of the lattice topology as well as the spatial 
dimensionality of the investigated system. It is shown that the diamond-like decorated planar lattices with a sufficiently 
high coordination number may exhibit a striking critical behaviour including reentrant phase transitions with two or 
three consecutive critical points. 



1 Introduction Quantum Heisenberg models on geo- 
metrically frustrated planar lattices have enjoyed a great 
interest during the past decades especially due to their ex- 
traordinary diverse ground-state behaviour, which is often 
a result of mutual interplay between geometric frustration 
and quantum fluctuations [ 1 , 2 3 1 . From this perspective, 
geometrically frustrated quantum systems represent an ex- 
cellent play ground for theoretical study of novel quantum 
many-body phenomena. Beside a rather complex ground 
state, which can be composed of several phases like the 
semi-classical Neel-like ordered phase, the quantum va- 
lence bond crystal phase or different disordered spin-liquid 
phases [1 1, quantum spin systems on two-dimensional ge- 
ometrically frustrated lattices furnish a deeper insight into 
the quantum order- from-disorder effect [fl~||2][3), the scalar 
or vector chirality [;4]Q, as well as, the non-zero residual 
entropy that characterizes the macroscopic degeneracy of 
the ground state [T) . 
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Among another interesting features, which also fre- 
quently attract the immense theoretical interest to the 
quantum Heisenberg model on geometrically frustrated 
planar lattices, belong an enhanced magnetocaloric effect 
emerging during the adiabatic demagnetisation |6]|2]|8] 
|9) and a presence of quantized magnetisation plateaux 
in low-temperature magnetisation curves Il2l l57l . It is 
worthwile to remark that both these outstanding physi- 
cal phenomena have also been observed in several real 
insulating magnetic materials representing possible ex- 
perimental realizations of two-dimensional geometrically 
frustrated Heisenberg models. Recent experimental stud- 
ies of the magnetocaloric effect in gadolinium galium 
garnet Gr3GasOi2 HI I 111 21 have revealed a potential ap- 
plicability of frustrated quantum spin systems as promis- 
ing refrigerant materials in the magnetic cooling tech- 
nique intended either for room temperature refrigera- 
tion or for satellite applications (see |6) and references 
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therein). The quantised magnetisation plateaux have been 
experimentally observed in the triangular lattice com- 
pounds CsFe(S0 4 ) 2 US, Cs 2 CuBr 4 fT4lp3tfT6tfT7]fT8| . 
and RbFe(Mo0 4 ) 2 p3l[l9p5||2T1|22|. the kagome lat- 
tice compound [Cu 3 (titmb) 2 (CH 3 COO) 6 ] -H 2 l23l . 
as well as the Shastry-Sutherland lattice compounds 
SrCu 2 (B0 3 ) 2 J24l|25l|26l|27l, RB 4 (R = Er,Tm) M 

Unfortunately, searching for the exact solution for the 
geometrically frustrated quantum Heisenberg models often 
fails due to a non-commutability between spin operators 
involved in their Hamiltonians. Owing to this fact, we have 
recently proposed a special class of geometrically frus- 
trated Ising-Heisenberg models on diamond-like decorated 
lattices l3^l33ll34li35ll36U37ll38ll . which can be examined 
within the framework of an exact analytical approach based 
on the generalised decoration-iteration transformation |39l 
|40PTI . These simplified quantum models overcome the 
afore-mentioned mathematical difficulty by introducing 
the Ising spins at nodal lattice sites and the Heisenberg 
dimers on interstitial decorating sites of the considered 
planar lattice. It is worth mentioning that the mixed-spin 
Ising-Heisenberg models with diamond-like decorations 
turn out to be very useful testing ground for elucidating 
many quantum properties of low-dimensional magnetic 
materials, in spite of the fact that the monomer-dimer 
interaction is considered as the Ising-type interaction. In- 
deed, these rather simple spin models shed light on diverse 
quantum features to appear in the ground state 132 3311341 
[35ll38l . the magnetisation process I34ll37ll38l . the thermo- 
dynamics IB4U37L as well as, the critical behaviour 1361 . 
Note that kinetically frustrated diamond chain models 
constituted by nodal Ising spins and mobile electrons de- 
localised over the interstitial decorating sites have been 
particularly examined as well 1421 14311441 . 

Motivated by these remarkable quantum features, be- 
ing not commonly observed in any semi-classical Ising 
system, we will investigate in the present work the mixed 
spin- 1/2 and spin-1 Ising-Heisenberg model on several 
decorated planar lattices consisting of inter-connected dia- 
monds. Beside the ground-state analysis, our main goal is 
to particularly examine critical behaviour of these quantum 
models. We will also endeavour to explore how a possible 
existence of the reentrant phase transitions depends on a 
particular choice of the lattice topology of the investigated 
spin system. 

The outline of this paper is as follows. In the next sec- 
tion, we will provide the detailed description of the Ising- 
Heisenberg model and then, some basic steps of the used 
decoration-iteration mapping method will be elucidated. 
Section [3] deals with the most interesting results for the 
ground-state and finite-temperature phase diagrams, which 
are supplemented by a brief discussion of temperature vari- 
ations of the total magnetisation and the specific heat. Fi- 
nally, some concluding remarks are drawn in Section|U 



2 Model and its exact solution Let us consider 
two-dimensional lattices composed of inter-connected dia- 
monds as is illustrated in Fig.[T]for honeycomb and square 
lattices. As shown in this figure, the diamond-shaped unit 
of considered lattices consists of two different kinds of spin 
sites schematically depicted as empty and full circles. To 
ensure an exact tractability of the investigated spin system, 
we will further suppose that the former spin sites are oc- 
cupied by the Ising spins a — 1/2 that interact with other 
spins through the interaction J\ containing just one (zth) 
spatial component of the spin operator, while the latter ones 
are occupied by the decorating Heisenberg spins 5 = 1 
that interact among themselves via the anisotropic XXZ 
coupling Jh (A) containing all three spatial components of 
spin operators. Note that the parameter A allows to con- 
trol the interaction Jh between the easy-axis (A < 1) and 
easy-plane (A > 1) type, as well as, to obtain the Ising 
model as a special limiting case when assuming A = 0. 
Under these assumptions, the total Hamiltonian of the in- 
vestigated model can be written as 

Nq/2 

7i = - J H ]T i A & § * + ) + 

2Nq 2Nq 

-J^SfK-D^iSf) 2 . (1) 

(i,n) 1=1 

Here, the first summation is carried out over all exchange 
interactions between pairs of the nearest-neigbouring 
Heisenberg spins, the second summation takes into account 
the interaction between the nearest-neigbouring Ising and 
Heisenberg spins and the last summation runs over all 
lattice sites occupied by the decorating Heisenberg spins. 
The spin variables S/ (7 = x, y, z) and er* denote spatial 
components of the spin- 1/2 and spin-1 operators located at 
ith and nth lattice position, respectively. The parameter D 
stands for the axial zero-field splitting (AZFS) parameter 
acting on the Heisenberg spins only |45 46 1. Finally, N 
represents the total number of the Ising spins and q is the 
coordination number of the lattice, which implies that the 
total number of all spins (lattice sites) is N to t = N(q+ 1). 

The crucial step of our calculation represents the eval- 
uation of the partition function of the system. In view 
of further manipulations, it is useful to rewrite the total 
Hamiltonian (Q~|) as a sum over the cluster Hamiltonians 

Tt = J2^2{ 2 ^fc> where each cluster Hamiltonian TCh in- 
volves all interaction terms that belong to just one (fcth) 
diamond unit (see lhs of Fig . O : 

H k = -MA0A + + S* kl S* k2 ] 

-MHi + faXHi + aU) - D[(S* kl ) 2 + (S z k2 n (2) 

Taking into account a validity of the commutation relation 
between cluster Hamiltonians of two different diamond 
units [Wfe,7ij] = 0, the partition function of the investi- 
gated mixed-spin Ising-Heisenberg model can be partially 
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(a) (b) 

Figure 1 The mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb [Fig. 02a)] and 
square [Fig. [Tib)] lattices. White circles represent the nodal Ising spins a = 1/2 and black ones denote the decorating 
Heisenberg spins S = 1. 




Figure 2 A schematic representation of the decoration- 
iteration transformation of the kth diamond-shaped unit of 
the investigated mixed-spin Ising-Heisenberg model (lhs of 
figure) towards the unit of the simple spin-1/2 Ising model 
on corresponding lattice with the nearest-neighbour inter- 
action R (rhs of figure). The parameters Jj, J^(A), and 
D stand for the exchange interaction between the nearest- 
neigbouring Heisenberg spins, the interaction between the 
nearest-neigbouring Ising and Heisenberg spins and the ax- 
ial zero-field splitting parameter, respectively. 

factorised into the following product of cluster partition 
functions Zk : 

Nq/2 Nq/2 

Z = E IT Tr fe exp(-/W fc ) = Y, IT Zfc < 3 > 

{ ai } k=l {at} k=l 

Here, (3 = l/(fcsT), fee is Boltzmann's constant and T is 
the absolute temperature. The summation ^{o- . } m Eq. (0 
runs over all possible spin configurations of the nodal Ising 
spins, while the symbol Tr^ stands for a trace over degrees 
of freedom of the fcth Heisenberg dimer. By performing 
an exact analytical diagonalisation of the cluster Hamilto- 
nian (0 in a particular Hilbert subspace corresponding to 
the Heisenberg spin pair from fcth diamond unit, the parti- 
tion function Zk will depend on two nodal Ising spins aki 
and <Jk2 only. Moreover, the resulting explicit form of Zk 
immediately implies a possibility of performing the gener- 



alised decoration- iteration mapping transformation I139II40I 

ETJ 

Z h = e 2 P D -? jH +2c 2t}D+l3jH cosh[2/?Ji(<r£ 1 + a z k2 )] 

+ 2c I 3d -P j k/ 2 cosh[/V(2D - J H ) 2 + 8(J H ^) 2 /2] 

+ Ae fiD cosh[/3Ji(cr^+cr^ 2 )]cosh(/3J H ^) 

= Ae^ 1 "' 2 . (4) 

From the physical point of view, the mapping transforma- 
tion removes all interaction parameters associated with 
the Heisenberg spins and replaces them by the effective in- 
teraction R between the remaining nodal Ising spins aki 
and <Tk2 (see Fig. [2] for better illustration). Of course, the 
transformation relation (0 must hold for any available spin 
configuration of the Ising spins. This "self-consistency" 
condition unambiguously determines so far not specified 
mapping parameters A and R, 

A = ^/WiW 2 , R = 2/r 1 (lnW 1 - lnW 2 ) . (5) 

It is clear that the functions W\ and W2 emerging in the 
mapping parameters (0 are two independent expressions 
to be obtained by substituting four possible configurations 
of the Ising spins into the cluster partition function ©: 

Wi = c 2130 - 13 ^ +2c 2i3D+i3Jh cosh(2/?Ji) 

+ 20^-^/2 CO sh[/V(2L> - Jh) 2 + S(JhA) 2 /2] 
+ 4e /3D cosh(/3Ji) cosh(/3J H zi), (6) 

W 2 = e 2 ^-^ j h + 2 c 2/5CI+ ' 3Jh 

+ 2e /3D -' 9jH / 2 cosh[/V(2£> - Jh) 2 + 8( J H Z\) 2 /2] 

+ 4c pD cosh(f3J H A). (7) 

At this stage, the straightforward substitution of the trans- 
formation relation (0 into the expression (0 yields the 
equality 

Z(T, J h J H , A, D) = A N ^ 2 Z lsing (T, R), (8) 

which establishes an exact mapping relationship between 
the partition function Z of the mixed spin-1/2 and spin-1 
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Ising-Heisenberg model on diamond-like decorated planar 
lattices and the partition function Zj S i ng of the simple spin- 
1/2 Ising model on corresponding undecorated lattices 
with the nearest-neighbour coupling R. In this respect, 
Eq. ([8]l formally completes our calculation for the parti- 
tion function of the spin- 1/2 and spin-1 Ising-Heisenberg 
model, since exact results for partition functions of many 
spin- 1/2 Ising planar lattices are known l40l . It is worth 
noticing that this mapping relation is universal and valid 
regardless of the lattice topology and space dimensional- 
ity of the investigated model system. Besides, it also al- 
lows a direct calculation of some physical quantities (such 
as Gibbs free energy, internal energy, magnetisation, en- 
tropy, susceptibility, or specific heat) by the use of stan- 
dard thermodynamical-statistical relations, which are im- 
portant for understanding of the magnetic behaviour of the 
investigated spin system. Other important physical quanti- 
ties, such as sub-lattice magnetisation mf and reduced 
per one Ising and Heisenberg spin, respectively, correla- 
tion functions cf ; z , cJJ, cgfj, and c-jf , as well as, quadrupo- 
lar moment q^, which cannot be obtained within afore- 
mentioned procedure, can be calculated by combining the 
relation (O with the exact mapping theorems developed by 
Barry et al. H47U48U49H and the generalized Callen-Suzuki 
spin identity I50H5TH521: 



mf = (a£i) = (a z kl ) = m , (9) 
mi={Sl 1 )=Am h{J l ), (10) 
c ii = fikivD = (°fcl<7fc2)o = £o, (11) 
tfi = = (l+4 £o )/ 2 (Ji) + (l-4e )/a(0),(12) 

ell = (S x kl S x k2 ) = (1+ 4e )/ 3 (Ji) + (1- 4e )/3(0),(13) 
cfh = (KJk) = (1/2 + 2£o)/i(Ji), (14) 



= ((S z kl ) 2 ) = (l+4 £o ).A(Ji) + (l-4e )/ 4 (0). (15) 

In above, the symbols (. . .) and (. . .)o represent stan- 
dard canonical averages performed over the ensemble 
defined by the mixed-spin Ising-Heisenberg model on 
the diamond-like decorated lattice and the spin- 1/2 Ising 
model on the corresponding lattice, respectively, mo labels 
the single-site magnetisation and £o stands for the nearest- 
neighbour correlation function of the corresponding Ising 
model. Finally, the functions fi(x) (i = 1-4), emerging in 
the set of Eqs. (II Oi l. dl21 — dl5l . are defined as follows: 

h(x) = [e^ D+/ " H sinh(2#r) 

+ ainh(j3x)coah(j3J K A)}/G(x), (16) 
f 2 (x) = [2c /3D+/3Jh cosh(2/3a;) - fJ> D -(> J * 

- (2D- J H )w- 1 c- /3jH/2 sinh(/3Lj/2) 

_ e - /3jH/2 cosh(/3w/2)]/[2G(a;)], (17) 
f 3 (x) = [cosh(/3x) sinh(/?J H Z\) 

+ 2J li Aoj- 1 e~ 0jH / 2 smh{f3uj/2)}/G(x), (18) 



U(x) = [2c /3D+/3jH cosh(2/3x) +c" /3jH / 2 cosh(/3w/2) 
+ {2D - Jh)w- 1 c- /3Jh / 2 sinh(/3w/2) + eP D ~^ 
+ 2 cosh(/3x) cosh(j3J K A)]/[2G(x)], (19) 

where G(x) = e 0D - 0jH + 2 e^ D+ ^ jH cosh(2/?x) + 
4cosh(/3x)cosh(/3J H Z\) + 2 e _/3jH/2 cosh(/3w/2) and 
UJ= y/(2D- J H ) 2 +8(J H /i) 2 . 

3 Results and discussion Before proceeding to 
a discussion of the most interesting results obtained for 
the mixed spin- 1/2 and spin-1 Ising-Heisenberg model on 
diamond-like decorated honeycomb and square lattices, let 
us make a few remarks on a validity of analytical results 
presented in the previous section. It is worth mentioning 
that all obtained results are universal as they hold regard- 
less of whether ferromagnetic or antiferromagnetic inter- 
action parameters Jj and Jh are assumed, as well as, in- 
dependently of the lattice topology or spatial dimension- 
ality of the investigated spin system. As proved, however, 
there are some fundamental differences between magnetic 
behaviour of models with distinct nature of the Heisen- 
berg interaction (see our preliminary reports I32il35ll36l ). 
Considering this fact, we will restrict ourselves here just 
to the case with the ferromagnetic Heisenberg interaction 
Jh > 0. To simplify our discussion, we will also assume 
the Ising interaction J\ to be ferromagnetic, because the 
sign change Jj — > —J\ brings just a trivial change in the 
local alignment of the nodal Ising spins with respect to 
their nearest Heisenberg neighbours. Indeed, it changes the 
parallel alignment between the Ising and Heisenberg spins 
along the z-axis to the antiferromagnetic one. Other par- 
ticular case of the system, in which both the interaction 
constants Jh and J\ are supposed to be antiferromagnetic 
(Jh < 0, Ji < 0), will be explored in the subsequent 
work [53). Bearing this in mind, we can set the Ising ex- 
change interaction J\ as the energy unit and introduce the 
following set of dimensionless parameters: t = h^T/ J\, 
d = D/J\, and a = Jh/Ji, as describing the dimension- 
less temperature, the relative strength of the AZFS param- 
eter, and the strength of the Heisenberg interaction normal- 
ized with respect to the Ising interaction, respectively. 

3.1 Ground-state properties Let us start our discus- 
sion with the analysis of the ground-state phase diagrams 
displayed in Fig. [3] As one can see from the phase diagram 
constructed in the A — lnaplane [Fig. Oa)], three different 
phases appear in total within the ground state of the sys- 
tem without AZFS parameter due to a mutual interplay be- 
tween the exchange interactions Jj, Jh, and the exchange 
anisotropy A. It is noteworthy that a presence of non-zero 
AZFS paramer does not lead to any new phase as illustrated 
on the ground-state phase diagram from Fig. Ob). Spin ar- 
rangements of three possible ground-state phases can be 
unambiguously characterised as follows: 
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Figure 3 Ground-state phase diagrams of the mixed spin- 
1/2 and spin-1 Ising-Heisenberg model with ferromag- 
netic Ising and Heisenberg exchange interactions (a) in 
the A — In a plane for the system without the AZFS pa- 
rameter and (b) in the A — d plane for the system with 
a = 1.0. Spin orders to emerge within the phase sectors 
FRO, QFO, and FRU are unambiguously determined by 
the wave-functions (I20li-(l22"li. 



N 



Nq/2 

iFRO)=ni+>«n i 1 ' 1 ! 



i=l 

mf = 0.5, m 



k=l 

he 



h — i., Cj; — 0.25, 
Ihh = 1> c hh = 1' c hh = 0) c ih = 0-5; 



(20) 



JV 



Nq/2 



iQFO)=ni+) i n ^(ii,o) fc +io,i) 

i=l k=l V 

mf = 0.5, mf = 0.5, cf z = 0.25, 

qli = 0.5, eg = 0, cH = 0.5, eff = 0.25; (21) 

N Nq/2 

ifru) = n i±>i n 2 hO 1 - - i)k + 1 - j > i)fe ) 

-V26_|0,0)J, 



fc=i 



mf = 0,mg = 0,cf z = 0, 

ffih = -cf h = 4/2, ell = 2J h Alj-\c% = 0, (22) 



where b± = y/l ± (2D — Jh)w _1 . In above, the first 
product is taken over all nodal Ising spins, while the second 
one runs over all Heisenberg dimers. Ket vectors |±) and 
| ± 1, 0) after relevant product symbols determine the spin 
states a z = ±1/2 and S z = ±1, of the Ising and Heisen- 
berg spins, respectively. Analytic expressions for the phase 
boundaries between these ground states read: 



|FRO) - |QFO) : d = a(A - 1) - 1, (23) 

A) 2 

(24) 



a (aA) 2 
|FRO)-|FRU): d=---l + ^—L- 
I Z(a + 1) 



| QFO) - | FRU) : 



a^ 2 
1+ 2, 



d- 



(25) 



which all meet at two triple points given by the condition 

Ti,2 = [At, dr] = 



1 + a VI + « VI + a 
~F j ~F 



(26) 



As one clearly sees from Eqs. ( f20l >, the phase FRO rep- 
resents the standard ferromagnetic phase, which can com- 
monly be observed in the pure Ising systems as well. Actu- 
ally, all the results indicate a perfect parallel alignment of 
the nearest-neighbouring Ising and Heisenberg spins in this 
phase. By contrast, other two phases QFO and FRU exhibit 
very unusual spin arrangements that cannot be observed 
in any semi-classical Ising systems. More specifically, all 
nodal Ising spins are perfectly aligned along the z-axis, 
while the decorating Heisenberg ones reside at a quan- 
tum superposition of spin states described by the symmet- 
ric wave function (| 1, 0) + |0, 1))/V2 in the QFO phase. 
Moreover, one may conclude from the location of the QFO 
phase in the phase diagram shown in Fig. [3ja) that this 
phase appears as a result of a mutual competition between 
the Ising interaction J\ favouring the ferromagnetic spin ar- 
rangement along the z-axis, and the easy-plane Heisenberg 
interaction Jn(A), which energetically favours the short- 
range ferromagnetic spin order in the xy plane. Of course, 
the position of QFO changes as the AZFS parameter d is 
switched on: for d > 0, one finds a lower range of val- 
ues of the exchange anisotropy A upon strengthening the 
AZFS parameter d corresponding to this phase, whereas 
this range at first extends and only then shrinks when d < 
[see Fig. EJb)]. Another phase FRU also appears as a re- 
sult of the mutual competition between the exchange in- 
teractions Ji and Jb_(A) with A > 1, respectively. As 
one can readily deduce from Eqs. ( l22l . this phase repre- 
sents a macroscopically degenerate monomer-dimer state, 
where all nodal Ising spins are frustrated due to antiferro- 
magnetic and/or 'non-magnetic' nature of the Heisenberg 
spin dimers. Actually, the Heisenberg spin pairs exhibit a 
quantum superposition of three spin states |0, 0), |1, —1), 
and | — 1, 1) in the FRU phase, whose relative probabilities 
depend on a mutual ratio between the interaction parame- 
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Figure 4 (a) The probability distribution for the entan- 
gled spin states of the Heisenberg spin pairs in the FRU 
phase as a function of the AZFS parameter d for A = 5.0 
and a = 1.0. (b) The corresponding zero-temperature de- 
pendencies of the quadrupolar moment ggg and correlation 
functions egg, egg. For clarity, the spin state with the major 
(minor) abundance is depicted in Fig.@ta) by full (broken) 
line. 

ters Jh, A, and D [see the eigenfunction in d22l l: 
P(\ ± 1, Tl» = -(l+ 2d ~ a ) , (27) 

P(|0, 0» = i ( 1 - ; 2d ~ a ) . (28) 

If 2p{\ ± 1, =Fl)) < p(|0, 0)), then the most probable spin 
state of the Heisenberg spin pairs is the 'non-magnetic' 
state 1 , 0) , while the antiferromagnetic | ± 1,^1) spin 
states become major microstates if 2p(| ± l,=pl}) > 
p(\0, 0)). According to this condition, the threshold value 
of the AZFS parameter dividing the parameter space of the 
FRU phase into a region with the dominant population of 
the 'non-magnetic' |0,0) spin state and a region with the 



prevailing population of the antiferromagnetic | ± 1, 
spin states is d* = 0.5a. 

More obvious insight into a probability distribution 
of both afore-mentioned kinds of microstates to emerge 
in the FRU phase may also be gained from Figs. Ua) 
and (b), which show the abundance probabilities for rel- 
evant spin states as functions of the AZFS parameter d 
for A = 5.0 and a = 1.0, as well as, the correspond- 
ing zero-temperature dependencies of the quadrupolar mo- 
ment q>gg and the correlation functions egg, egg, respec- 
tively. Note that the curves corresponding to the probabil- 
ities p(\ ± 1, Tl)) ar, d p(|0, 0)) are for clarity depicted in 
Fig- Eta) as full (broken) lines when the spin configuration 
| ± 1, =pl) and |0, 0) of the Heisenberg spin pairs has the 
major (minor) abundance in the FRU phase. As one can 
see from this figure, the probability to find the Heisenberg 
spin pairs in the 'non-magnetic' 1 0, 0) spin state gradually 
decreases upon increase of the AZFS parameter, while the 
probability of finding two antiferromagnetic | ± 1, =pl) spin 
states gradually increases with increasing the AZFS pa- 
rameter d. As a result, the abundance of these two kinds 
of spin states changes as soon as the threshold value of the 
AZFS parameter d* is attained. Above this boundary value, 
the population of the antiferromagnetic | ± 1, =pl) ('non- 
magnetic' |0,0)) spin states becomes major (minor) and 
the probability of its finding further increases (decreases) 
as d increases. 

As an independent check of this scenario, we have 
plotted the zero-temperature dependencies of the pair cor- 
relation function egg and the quadrupolar moment ggg 
in Fig. HJb). Obviously, both these physical quantities as 
functions of the parameter d exhibit monotonous depen- 
dencies: the correlation function egg monotonically de- 
creases, while the quadrupolar moment ggg monotonically 
increases upon increase of d. By contrast, the correlation 
function egg exhibits a much more complex variation: 
it first increases to reach its maximum at d* and then 
decreases with the increasing d. As expected, this non- 
monotonic behaviour of egg relates to a different trend of 
the short-range order of the Heisenberg spins in the region 
d < d* with the most populated 'non-magnetic' |0, 0) 
spin state and the region d > d* with the most populated 
antiferromagnetic | ± 1,^1) spin states. In the former re- 
gion, the increase in the AZFS parameter d reinforces the 
influence of the easy-plane exchange anisotropy A > 1. In 
the consequence of that, the short-range antiferromagnetic 
order along the z-axis is accompanied with the short-range 
ferromagnetic order in the xy plane. The precisely opposite 
situation emerges in the latter region: for any d < d* the 
exchange anisotropy A > 1 is restrained by the parameter 
d, which leads to the destruction of the short-range ferro- 
magnetic order of the Heisenberg spins in the xy plane on 
behalf of short-range antiferromagnetic spin order along 
the z-axis whenever d increases. 

3.2 Finite-temperature behaviour To provide a 
deeper insight into the finite-temperature behaviour of 
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Figure 5 The critical temperature of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated 
honeycomb [Figs.|3a), (b)] and square [Figs.|3c), (d)] lattices at the fixed interaction ratio a = 1.0 as a function of the 
exchange anisotropy A for several values of the AZFS parameter d and as a function of the AZFS parameter d for several 
values of the exchange anisotropy A. 



the investigated model system, let us turn our attention to 
the phase diagrams displayed in Fig. [5] This figure shows 
the critical temperature of the mixed spin-1/2 and spin-1 
Ising-Heisenberg model on diamond-like decorated hon- 
eycomb and square lattices as a function of the exchange 
anisotropy A and the AZFS parameter d for the fixed inter- 
action ratio a = 1.0. As it can be clearly seen, the model 
under investigation displays a rather complex critical be- 
haviour that very sensitively depends on a strength of both 
parameters A and d, as well as, the topology (coordina- 
tion number) of the investigated lattice. Apart from the ex- 
pected decrease of t c with the increasing A and/or decreas- 
ing d, one may also find non-monotonic dependencies of t c 
to be closely related to the FRO — > QFO phase transition 
[see for instance the curves labeled as d = —1.4, —1.45 
in Fig. |3 a) and the curves A = 3.3, 3.5 in Fig. HJb)]. 
Moreover, several interesting regions with reentrant phase 
transitions can also be observed in phase diagrams of the 
diamond-like decorated square lattice. More specifically, 
the square lattice shows the reentrant behaviour with ei- 
ther two or three consecutive critical points [see the curves 
d = -1.45,-1.47 in Fig. 0c) and the curve A = 3.5 



in Fig. |5jd)] . According to the ground-state analysis, this 
lattice exhibits two or three consecutive phase transitions 
of second order only if the anisotropy parameter A and 
the AZFS parameter d are selected sufficiently close to the 
phase transitions FRO-FRU and QFO-FRU. In general, 
the origin of this non-trivial phenomenon lies in a mutual 
competition between the easy-plane (D < 0, A > 1) and 
easy-axis (D > 0, A < 1) interactions, while the latter 
one is also supported by the Ising interaction J\. Evidently, 
the reentrance occurs either if the AZFS parameter d takes 
the negative values and the exchange anisotropy A is of 
an easy-axis type or if the reverse case is considered [see 
the curves d = —1.45, —1.47 and also the curve A = 3.5 
in Figs.|5jc) and (d), respectively]. By comparing the for- 
merly published results for the diamond-like decorated 
triangular lattice (see Ref. [36]) with the results shown 
in Fig. [5] one may conclude that the parameter region 
with reentrant phase transions enlarges with the increasing 
coordination number of the lattice. Indeed, the mixed spin- 
1/2 and spin-1 Ising-Heisenberg model on the decorated 
triangular lattice exhibits two or three critical points also 
for D < and A > 1. Besides, this lattice as the only one 
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is characterised by very interesting dependencies of t c that 
have a shape of semi-closed and closed loops (see Fig. 2 
in Ref. 11361 ). In these parts of phase diagrams, the system 



§ 0.4 











a = 


1.0 










d = 


-1.47- 




s/1 = 0.0 

. o — — 




<? = 


4 




1/ 

/ 2.25 




j_ K> : 


0.4 










1 O : 


o 









°0.0 0.1 0.2 0.3 0.4 0.5 



t 

Figure 6 Temperature dependencies of the total magneti- 
sation of the decorated square lattice with the interaction 
ratio a = 1.0 and the AZFS parameter d = —1.47 for 
several values of the exchange anisotropy A. 

starts from the disordered ground state before entering the 
spontaneously ordered phase (either FRO or QFO) at lower 
critical temperature t c i, which very quickly disappears due 
to strong thermal fluctuations at upper critical temperature 

tc2- 

The results displayed in Fig.[5]can be convincingly ev- 
idenced by thermal dependencies of the total magnetisa- 
tion, which represents the order parameter for both ferro- 
magnetic phases FRO and QFO. For this purpose, some 
temperature variations of the total magnetisation m = 
(raf +4m^)/5 of the diamond-like decorated square lattice 
are displayed in Fig.|6]for the interaction ratio a = 1.0, the 
AZFS parameter d = —1.47 and several values of the ex- 
change anisotropy A. Evidently, the m(t) curve exhibits 
reentrant transitions with two or three consecutive criti- 
cal points for A = 0.4 and 0.25, respectively, in accor- 
dance with the finite-temperature phase diagram displayed 
in Fig. SJc). Note that both these values of the exchange 
anisotropy relate to the regions rather close to the phase 
transition between the ferromagnetically ordered phase and 
the disordered FRU phase. On the other hand, when one is 
considering the exchange anisotropies that are sufficiently 
far from the phase transitions FRO-FRU and QFO-FRU, 
then the plotted magnetisation curves m(t) exhibit R-type 
dependencies with a single critical temperature (see e.g. the 
curves A = 0.0 and 1.0). According to the ground-state 
analysis, the initial values m = 0.9 and 0.5 correspond to 
the FRO and QFO phases, respectively [for clarity see also 
the set of Eqs. ® and dUll. 

Finally, let us close our discussion by exploring tem- 
perature dependencies of the specific heat. For illustation, 
some typical thermal variations of the specific heat of the 
Ising-Heisenberg model on the diamond-like decorated 



square lattice are plotted in Fig. [7] To enable a direct com- 
parison, we have chosen the values of the interaction ratio 
a, the AZFS parameter d and the exchange anisotropy 
parameter A so as to match all temperature dependencies 
of the total magnetisation plotted in Fig. [6] In this manner, 
the displayed set of specific heat curves reflects a com- 
prehensive picture of the finite-temperature behaviour of 
the investigated spin system. If the exchange anisotropy 
parameter A is chosen to be sufficiently far from the phase 
transitions FRO-FRU and QFO-FRU, then temperature 
dependencies of the specific heat with just one logarithmic 
singularity are observable at the continuous (second-order) 
phase transition between the spontaneously ordered and 
disordered phases [see Figs. [3 a) and (d)]. It should be also 
mentioned that these thermal dependencies exhibit yet an- 
other two or three round Schottky-type maxima in addition 
to a pronounced logarithmic divergence. On the other hand, 
remarkable specific heat curves with two or three logarith- 
mic singularities can also be detected when the value of A 
is taken from the region close to the phase transitions FRO- 
FRU and QFO-FRU [see Figs. 0b) and (c)]. Among other 
matters, the observed logarithmic singularities provide an 
independent check of the reentrant behaviour previously 
discussed by finite-temperature phase diagrams, as well 
as, temperature dependencies of the total magnetisation. 
Besides, one or two local maxima can also be detected in 
low-temperature tails of all displayed specific heat curves 
in addition to a broad Schottky-type maximum to emerge 
at high temperatures. It is quite reasonable to expect that 
the origin of observed low-temperature maxima lies in 
strong thermal excitations to a spin configuration rather 
close in energy to the ground state. 

4 Concluding remarks In the present paper, the 
ground-state and critical behaviour of the mixed spin- 
1/2 and spin-1 Ising-Heisenberg model on diamond-like 
decorated honeycomb and square lattices has been investi- 
gated within the framework of the generalised decoration- 
iteration mapping transformation. Using this rigorous pro- 
cedure, the exact solution for the investigated mixed-spin 
model has been obtained by establishing a precise mapping 
equivalence with the spin- 1/2 Ising model on the corre- 
sponding undecorated planar lattice with the known exact 
solution. 

The main emphasis of this work was to bring a deeper 
insight into the effect of quantum fluctuations and geo- 
metric frustration generated by the mutual competition 
between the Heisenberg interaction, Ising interaction and 
AZFS parameter on magnetic properties of the planar 
mixed-spin model with diamond-like decorations at zero 
as well as non-zero temperatures. Our further goal was 
to shed light on how the critical behaviour of this sys- 
tem depends on its lattice topology. The obtained results 
clearly demonstrate that an interplay between the Ising in- 
teraction favouring the spin arrangement along the z-axis 
and the easy-plane XXZ Heisenberg interaction favouring 
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the short-range ferromagnetic order in the xy plane gives 
rise to two interesting ground states (QFO and FRU) with 
entangled states of the Heisenberg spins. Moreover, the 
appearance of one of those phases (namely, the appearance 
of the FRU phase) verified the known fact that the compe- 
tition between the ferromagnetic Ising and ferromagnetic 
easy-plane Heisenberg interactions leads to occurrence 
of the short-range antiferromagnetic spin order along the 
z-axis, which is accompanied by the short-range ferro- 
magnetic spin order in the xy plane [32,54,55]. Hence, 
one may conclude that this behaviour represents a general 
feature of all quantum models with a mixture of easy-axis 
and easy-plane bonds, since it can be found regardless of 
the lattice topology as well as the spatial dimensionality of 
the investigated model system. 

The most interesting result to emerge from the present 
study consists in the exact evidence of the existence of an 
reentrant phase transitions with either two or three consec- 
utive critical points. The existence of this non-trivial phe- 
nomenon has also been evidenced by temperature depen- 
dencies of the total magnetisation, as well as, remarkable 
thermal variations of the specific heat exhibiting two or 
three logarithmic singularities. As it has been proved, the 
higher is the coordination number of the lattice, the more 
pronounced the reentrance is and more diverse critical be- 
haviour of the system may be found. 

Finally, it is worthwhile to remark that even though our 
theoretical investigation of the mixed spin-1/2 and spin- 
1 Ising-Heisenberg model on diamond-like decorated pla- 
nar lattices was mainly aimed at providing a deeper insight 
into cooperative and quantum features of this exactly solv- 
able model, we hope that our results might stimulate re- 
searches on possible experimental realizations of this in- 
teresting quantum spin model and confirm our theoreti- 
cal predictions. From this perspective, the most promising 
approach in experimental realizations of our model sys- 
tem may represent a targeted design of cyano-based poly- 
meric coordination compounds or their isostructural ana- 
logues. For instance, a series of bimetallic polymeric coor- 
dination compounds Cu(L) 3 [Fe(CN) 6 ](C104) 2 -nH 2 0, L = 
N-(3-aminopropyl)-l,3-propanediamine 11561571 or N-(2- 
aminoethyl)-l,3-propanediamine [58], with the diamond- 
like decorated honeycomb network structure may repre- 
sent useful starting point for this rational synthesis. In 
this series, the divalent Cu 11 and Fe 11 metal atoms re- 
side decorating and nodal sites of the diamond-like dec- 
orated honeycomb lattice, respectively [Fig. HI a)]. Unfor- 
tunately, the divalent Fe 11 atoms are due to a strong ligand 
field of the cyano group in the diamagnetic low-spin state 
with 5 = 0. Similarly, the bimetallic polymeric coordi- 
nation compounds [W(CN) 4 Fe (H 2 0) 22 ]nH 2 J59) and 
[W(CN) 4 Co(H 2 0) 22 ] nH 2 gO) with the diamond-like 
decorated square network structure [Fig. QIb)] have also 
been reported quite recently, but here, the tetravalent W IV 
metal atoms residing nodal sites of the diamond-like dec- 



orated square lattice are due to a strong ligand field of the 
cyano group in the diamagnetic state 5 = 0. 
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